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Purpose:  This  Coastal  Engineering  Technical  Note  provides  guidance  on  using  movable-bed 
physical  models  to  predict  erosion  and  deposition  caused  by  currents  in  tidal  inlet  channels.  A 
scaling  relationship  based  on  equilibrium  scour  depth  allows  observed  model  scour  depths  to  be 
quantitatively  scaled  to  full-scale  dimensions.  Appropriate  situations  are  listed  for  which  this 
modeling  guidance  can  be  applied. 

Background:  Historically  movable-bed  physical  models  of  jettied  and  jettied  inlet  systems 
have  seen  limited  use  due  to  uncertainty  about  similitude  relationships  for  scaling  model  results 
to  full  size  (prototype  scale).  Consequently,  most  movable-bed  inlet  model  results  were  viewed 
as  qualitative  indicators  of  general  inlet  behavior  and  evolution,  and  little  attempt  has  been  made 
to  use  movable-bed  models  to  quantify  depth  of  scour  or  amount  of  deposition  under  given  flow 
conditions. 

Area  constraints  in  laboratory  facilities  generally  force  the  prototype-to-model  geometric  length 
scale  (Nl)  of  a  typical  inlet  to  be  on  the  order  of  75-to-l  or  greater  (NL  >  75).  For  a  model  in 
which  the  flow  kinematics  conform  to  the  Froude  similitude  relationship,  it  is  not  practical  to 
reduce  the  model  sediment  size  by  the  geometric  length  scale  and  still  have  noncohesive 
sediment  in  the  model.  Instead,  noncohesive  sediment  is  used  in  the  model  that  has  a  grain-size 
diameter  relatively  larger  than  the  diameter  required  by  Froude  scaling.  The  consequence  is  that 
model  currents  need  to  be  relatively  faster  than  currents  specified  by  Froude  scaling  to  move  the 
sediment  and  scour  the  inlet  channel  to  a  depth  similar  to  what  would  occur  in  nature.  Therefore, 
it  is  necessary  to  “distort”  the  prototype-to-model  Froude  velocity  scale  to  achieve  similarity  in 
scour  patterns  between  the  model  and  prototype.  Unfortunately,  distorting  the  velocity  scale 
limits  the  movable-bed  model  to  situations  where  inlet  channel  scour  is  caused  by  either  currents 
or  waves,  but  not  by  both. 

This  note  focuses  on  movable-bed  modeling  of  inlet  channel  scour  problems  stemming  from 
bedload  sediment  transport  caused  by  tidal  currents.  Such  scour  is  most  common  in  the  throat  of 
inlet  systems  where  waves  have  only  a  secondary  effect  because  of  reduced  wave  height  or  fairly 
deep  channel  depths.  Situations  where  waves  are  thought  to  be  a  dominant  scour  mechanism  are 
not  correctly  simulated  by  the  relationship  given  in  this  note. 

Scaling  Relationship:  An  appropriate  scaling  relationship  for  movable-bed  modeling  of 
channel  evolution  produced  by  the  tidal  current  is  derived  in  Appendix  A  of  this  note.  The 
scaling  relation  was  derived  by  assuring  that  an  expression  for  equilibrium  scour  depth  at  tidal 
inlets  is  the  same  in  the  model  as  in  the  prototype.  The  resulting  scaling  relationship  is  given  by 
the  expression: 
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in  which  the  scale  ratios  are  defined  as: 

N-  =  velocity  scale  [=Vp/vmVp  Nm\ 

N(Ss-i)  =  sediment  immersed  specific  gravity  scale  [=(Ss  -  l)p/(Ss  -  7)m] 

Nde  =  sand  grain  size  scale  [=(de)vl{de)m] 

Nl  =  geometric  length  scale  [=L//L/„J 

where  the  subscripts  p  and  m  represent  “prototype”  and  “model,”  respectively,  and  the  variables 
are  defined  as: 


V  =  depth-averaged  velocity 

Ss  =  sediment  specific  gravity  [=  ps  / 'pw ] 

ps  =  mass  density  of  sediment 

pw  =  mass  density  of  water 

de  =  sand  median  grain  size 

L  =  length 

The  N(Ss-i)  scale  ratio  allows  use  of  model  sediment  having  different  density  than  the  prototype. 

However,  for  practical  reasons  quartz  sand  is  typically  used  in  models  so  N(ss-i)  =  1,  reducing 

Equation  1  to  a  simpler  scaling  relationship  given  by: 

N-=N^nT  (2) 


The  Nde  scale  ratio  in  Equation  2  compensates  for  model  sediment  that  is  relatively  larger  than  if 
it  had  been  geometrically  scaled  according  to  the  length  scale,  Nl,  as  required  by  strict  geometric 
scaling.  Consequently,  the  velocity  scale  in  Equations  1  and  2  is  a  “distorted”  velocity  scale.  In 
other  words,  the  sediment  grains  in  the  model  are  larger  than  they  should  be,  so  the  model 
currents  need  to  be  faster  than  Froude-scaled  currents  to  achieve  the  same  equilibrium  depth  of 
scour.  For  the  extremely  rare  case  when  the  prototype  sediment  can  be  scaled  geometrically, 
then  Nde  =  Nl,  and  Equation  2  reduces  to  the  Froude  velocity  scale,  given  by 

N-  =  <  (3) 

Scaling  Validation:  The  scaling  relationship  given  by  Equation  2  is  validated  by 
demonstrating  that  results  from  small-scale  laboratory  inlets  at  (or  approaching)  equilibrium 
follow  the  same  tidal-prism  relationship  as  real  inlets  when  the  model  results  are  scaled  using 
Equation  2. 

O’Brien  (1931,  1969)  established  the  relationship  between  inlet  equilibrium  cross-sectional  area 
and  tidal  prism.  O'Brien’s  original  work  was  extended  and  refined  by  Jarrett  (1976).  The  solid 
line  on  Figure  1  is  Jarrett’ s  Area  -  Prism  relationship  for  all  inlets.  The  adjacent  dashed  lines  are 
the  95  percent  confidence  limits.  The  cluster  of  data  points  in  the  lower  left  comer  of  Figure  1 
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Figure  1 .  Model  inlet  results  scaled  to  prototype  (To  obtain  meters,  multiply  feet  by  0.3048) 


are  results  from  18  small-scale  movable-bed  model  experiments.  The  solid  circles  are  11 
experiments  conducted  by  Mayor-Mora  (1977),  and  the  solid  triangles  are  7  experiments 
conducted  as  part  of  the  Coastal  Inlets  Research  Program  (CIRP).  These  experiment  results  fall 
well  below  the  tidal  prism  relationship  indicating  that  the  inlet  cross-sectional  area  in  the  model 
was  an  order  of  magnitude  larger  than  what  might  be  expected  for  the  model  tidal  prism. 

Model  results  were  scaled  to  prototype  size  using  both  conventional  Froude  scaling  and  the 
scaling  relationship  given  in  this  technical  note.  For  both  methods  the  tidal  prism  can  be 
represented  as: 

p  _  Qmax  T  _  A  Vfnc/x  T 
n  k 

where 

P  =  tidal  prism 

Qmax  =  maximum  flow  discharge 
T  =  tidal  period 
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A  =  inlet  throat  cross-sectional  area 
Vmax  =  average  maximum  velocity  at  inlet  throat 

The  tidal  prism  scale  ratio  is  given  by  the  prototype-to-model  ratio  of  Equation  4,  i.e.,: 

N p  =  Nv  NT  N A  (5) 

where  Ny,  Nt,  and  NA  are  the  velocity  scale,  tidal  period  time  scale,  and  cross-sectional  area 
scale,  respectively. 

Froude  Scaling.  Noting  for  Froude-scaled  models  Ny  =  Nt,  Na  =  (NLf,  and  Nl  =  (Nt)2,  the 
prism  scale  can  be  expressed  in  terms  of  just  the  time  scale.  Thus,  the  model  results  can  be 
scaled  to  prototype  area  and  prism  using  the  Froude-scale  equations: 

Ap  =  Am{NT)4  and  Pp  =  Pm(NT)6  (6) 

where  once  again  the  subscripts  p  and  m  represent  prototype  and  model,  respectively.  Model 
tidal  periods  were  used  to  establish  the  time  scale,  and  the  corresponding  length  scales  were  then 
fixed  by  the  Froude  scale.  (Resulting  prototype  length  scales  ranged  between  50  <  NL<  1,388.) 
The  Froude-scaled  experiment  results  are  shown  as  x-symbols  (Mayor-Mora  1977)  and  plus- 
signs  (CIRP  data)  on  Figure  1.  All  but  one  value  are  outside  Jarrett’s  95  percent  confidence 
limit,  with  cross-sectional  areas  smaller  then  would  be  expected  for  the  scaled-up  tidal  prism. 
This  is  logical  because  the  model  sediment  diameter  was  much  larger  than  required  for  Froude 
scaling,  so  it  was  not  as  easily  scoured  from  the  model  inlet. 

New  Scaling  Relationship.  Scaling  the  small-scale  model  results  to  prototype  using  the 
scale  relationship  given  in  this  technical  note  uses  the  same  tidal  prism  scale  given  by 
Equation  5.  But  instead  of  Froude  velocity  scaling,  Equation  2  is  substituted  for  Ny  in 
Equation  5,  yielding 

np  =  Nd*  nt  nl'S  (7) 

The  relationship  between  the  length  and  time  scales  is  found  from  the  identity  Ny  =  Nl/Nt  or: 

Nl  =  Nf  Nj1  (8) 

u,e 

where  Ny  is  given  by  Equation  2.  Substituting  Equation  8  into  Equation  7  results  in  the 
following  relationships  for  scaling  area  and  prism  to  prototype  scale: 

Ap  =  Am{NL)2  =  Am(Ndefn  (NTfn  and  Pp  =  Pm  (Nif1  {Nt)24P  (9) 

The  model  sediment  grain  size  was  0.34  mm  in  Mayor-Mora’ s  (1977)  experiments,  and  the  CIRP 
experiments  used  grain  sizes  of  0.13  mm  and  0.26  mm.  Because  0.34  mm  is  typical  for  real  tidal 
inlets,  this  value  was  chosen  as  the  prototype  grain  size.  Therefore,  Nde  was  unity  for  Mayor- 

4 


ERDC/CHL  CETN-IV-26 
June  2000 

Mora’s  data  and  either  1.3  or  2.6  for  the  CIRP  data.  Using  the  same  time  scales  as  before,  CIRP 
and  Mayor-Mora’s  (1977)  model  results  were  scaled  to  prototype  and  plotted  as  hollow  circles 
and  triangles  in  Figure  1.  (Resulting  prototype  length  scales  ranged  between  10  <  Nl  <  70.) 
Although  not  perfect,  most  points  fall  within  the  95  percent  confidence  band  associated  with  the 
tidal-prism  relationship,  indicating  reasonable  agreement  with  real  inlets.  Also  note  that 
sediment  grain  size  was  included  in  the  scaling,  although  its  influence  is  relatively  minor. 

Appropriate  Movable-Bed  Modeling  Applications:  The  most  critical  aspect  of  any 
modeling  technology  is  understanding  which  real-world  situations  are  appropriate  candidates  for 
modeling.  This  is  particularly  true  for  movable-bed  physical  models  because  of  potential  “scale 
effects”  related  to  the  model  sediment. 

The  scaling  relationship  given  in  this  note  reduces  sediment  scale  effects  by  attempting  to  assure 
that  the  balance  between  the  boundary  layer  shear  stress  acting  on  the  bottom  and  the  critical 
shear  stress  of  the  bed  material  is  preserved  in  the  model.  This  “shear  stress  balance”  restricts 
the  modeling  technology  to  portions  of  real  inlets  with  the  following  characteristics: 

•  Bottom  scour  is  primarily  due  to  the  tidal  current 

•  Concurrent  wave  action  is  small  and  does  not  contribute  significantly  to  sediment 
transport 

•  Sediment  is  transported  in  bed-load  mode  (approaching  equilibrium) 

•  Sediment  is  noncohesive  with  only  minor  cohesive  components 

Typical  inlet  regions  and  processes  that  could  be  modeled  using  the  scaling  guidance  presented 
in  this  technical  note  are  listed  as  follows: 

Inlet  Throat  Section.  All  or  portions  of  a  structured  or  unstructured  inlet  throat 
section  could  be  modeled.  However,  attempting  to  model  entire  throat  sections  for 
large  inlets  is  not  practical  without  introducing  geometric  distortion.  (Geometric 
distortion  is  possible  using  the  proposed  scaling  guidance,  but  several  other  factors 
must  be  considered  as  well.) 

Effects  of  Flow  Jets.  Jet  flow  features  created  by  inlet  jetty  planform  geometry 
will  be  reproduced  by  the  model,  and  the  model  bed  evolution  will  respond 
accordingly.  However,  the  relatively  faster  model  velocity  of  the  jet  will  enhance 
flow  entrainment  at  the  jet  boundary,  and  this  effect  should  be  analyzed. 

Effects  of  Structures.  Localized  bed  evolution  adjacent  to  the  channel  side  of 
jetty  structures  will  be  correctly  simulated,  as  will  scour  at  the  tip  of  training 
structures.  Scour  at  free-standing  bridge  piers  will  probably  not  be  in  similitude 
because  the  dominant  horseshoe  vortex  causing  scour  at  vertical  piers  and  piles  is  a 
different  mechanism  than  the  shear-stress  balance  assumed  in  this  scaling 
relationship  derivation. 
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Inlet  Modifications.  Effects  to  channel  depth  brought  about  by  jetty  modification 
and  extension,  construction  of  additional  structures,  and  navigation  channel 
realignment  can  be  simulated.  However,  regions  of  deposition  will  require  a  source 
of  infilling  sediment. 

Modeling  Caveats:  Two  caveats  apply  to  movable-bed  modeling  of  the  situations  described 
above.  First,  this  modeling  technology  pertains  only  to  the  equilibrium  condition  where  the 
movable  bed  no  longer  evolves.  Sediment  transport  rates  and  the  time  required  to  reach 
equilibrium  are  not  correctly  simulated  in  the  movable-bed  model.  In  other  words,  the  model 
will  provide  reasonable  predictions  of  the  final  equilibrium  bed  configuration,  but  no  information 
is  gained  regarding  how  long  this  evolution  will  take  in  the  real  world.  This  is  actually  beneficial 
because  it  precludes  having  to  simulate  tidal  cycles  in  the  model.  At  real  inlets,  the  final 
equilibrium  condition  is  maintained  for  the  most  part  by  the  peak  ebb  or  flood  flow.  During  the 
rest  of  the  tide  stage,  currents  are  reduced  and  sediment  motion  is  greatly  reduced.  Therefore,  in 
the  movable-bed  model  we  only  need  to  simulate  the  maximum  flow  condition  long  enough  to 
achieve  an  equilibrium. 

The  second  caution  applies  to  deposition  areas  in  the  movable-bed  model.  Sediment  deposition 
will  occur  in  areas  where  the  local  flow  velocity  falls  below  the  sediment  incipient  motion 
threshold.  Sediment  moved  into  this  region  is  deposited.  The  scaling  premise  of  the  movable- 
bed  model  scaling  relationship  covers  this  situation,  but  there  must  be  available  sediment  to  have 
deposition.  At  real  inlets,  sediment  sometimes  comes  from  regions  of  the  inlet  that  are  actively 
scouring  because  of  structural  modification  or  change  in  flow  condition.  However,  most 
sediment  moving  through  an  inlet  is  being  continually  introduced  via  the  longshore  drift.  Some 
sediment  swept  into  the  inlet  deposits  on  the  flood  shoal,  some  deposits  in  the  dredged 
navigation  channel,  and  some  moves  out  of  the  inlet  to  the  ebb  shoal.  Therefore,  if  a  sediment 
deposition  problem  such  as  channel  infilling  is  to  be  studied  in  a  movable-bed  model,  it  is 
necessary  to  introduce  the  correct  quantity  of  sediment  into  the  model  at  the  appropriate 
locations. 

The  best  modeling  practice  is  to  validate  the  movable-bed  model  by  reproducing  the  existing 
condition  before  exploring  design  alternatives.  This  can  only  be  done  for  regions  of  inlets  that 
are  close  to  equilibrium.  Provided  the  movable-bed  model  bed  evolution  is  in  reasonable 
agreement  with  the  prototype,  then  the  effects  of  inlet  modifications  can  be  studied  with  greater 
confidence  that  the  model  results  are  good  predictors  of  the  expected  behavior  of  the  real  inlet. 
Hughes  and  Schwichtenberg  (1998,  1999)  successfully  validated  a  movable-bed  model  of  scour 
that  occurred  on  the  lee  side  of  the  Ventura  Harbor  breakwater.  Once  validated,  the  model  was 
used  to  design  suitable  scour  protection  for  the  breakwater  toe. 

Example  Application:  Determine  the  parameters  for  a  movable-bed  model  of  the  entrance 
channel  at  Shinnecock  Inlet,  New  York.  The  size  of  the  model  facility  requires  the  prototype- to- 
model  scale  ratio  to  be  Nl  =  75,  and  the  model  sediment  is  quartz  sand  with  a  median  grain  size 
of  (de) m  =  0.13  mm. 
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Prototype  Parameters: 

Inlet  width:  Wp  =  245  m 
Maximum  depth:  {hmax)v  ~  1 4  m 
Median  grain  size:  (de)p  =  0.61  mm  (quartz) 

Maximum  velocity:  V p  ~1.8  m/s  (depth- averaged) 

Maximum  discharge:  (Qmax) p  ~  2,400  m3/s 

The  inlet  width  and  maximum  depth  are  geometric  parameters  that  scale  directly  with  the  length 
scale,  e.g.,: 


W 

vyp 

W 

r  v  m 


245 


m 


=  N  t  =75  or  Wm  =— * ^  =  :i^:  =  3.3m(l0.7ft) 
L  m  75  75  - V  ’ 


The  grain-size  diameter  scale  ratio  is  given  as: 


N w  = 


M, 


0.61mm 
0.13  mm 


=  4.7 


which  is  needed  to  determine  the  required  velocity  scale  from  Equation  2,  i.e.,: 


N-  =  Nj*  N1^  =  (4.7) 3/8  (75) 1/8  =  107 


The  maximum  model  velocity  comes  directly  from  the  velocity  scale  as: 

£^  =  3.07  or  V m  «  «  1-8  m/S  «  0.59  m/s  (1.9  ft/s) 

Vm  3.07  3.07  - 1  ’ 

Recognizing  that  the  discharge  scale  is  simply  the  velocity  scale  multiplied  by  the  area  scale, 
i.e.,: 

Nq  =  Ny  Na  =  Ny  nI  =  3.07  (75)2  =  17,270 
the  maximum  model  discharge  is  obtained  as: 


— =  17,270  or  - 

Qm  17,270  17,270 


2-400m/s  0.14m3/s(4.9ft3/s) 


In  summary,  the  movable-bed  model  of  Shinnecock  Inlet  should  have  the  following 
characteristics.  These  are  reasonable  parameters  well  within  laboratory  capabilities. 
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Model  Parameters 


Inlet  width:  Wm 
Maximum  depth  (hmax)in 
Median  grain  size:  (de)in 
Maximum  velocity:  Vm 
Maximum  discharge  (Qmax)m 


=  3.3  m 
*0.19  m 

=  0.13  mm  (quartz) 

*0.59  m/s  (depth  averaged) 
*0.14  m3/s 


NOTE:  If  model  velocities  are  scaled  according  to  Froude  scale,  the  required  model  velocity  would  be: 

-1^=1^ 
m  a/75 


which  is  beneath  the  critical  velocity  necessary  to  initiate  movement  of  0.13  mm  quartz  sand. 


Additional  Information.  Questions  about  this  note  can  be  addressed  to  Ms.  Jackie  Pettway, 
Jackie.S.Pettway@usace.amiv.mil  This  note  was  produced  under  the  Coastal  Inlets  Research 
Program  (CIRP).  For  further  information  on  the  CIRP,  please  contact  the  CIRP  Technical 
Leader,  Dr.  Julie  Rosati,  Julie.Rosati@usace.anTiv.mil .  This  note  should  be  cited  as  follows: 


Hughes,  S.  A.  (2000).  "Simulation  of  current-induced  scour  in  movable-bed  inlet  models," 
ERDC/CHL  CETN-IV-26,  U.S.  Army  Engineer  Research  and  Development  Center,  Vicksburg, 

MS.  (httpf/chl. wes, army, mil/librarv/vublications/cetn) 
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Appendix  I:  Derivation  of  Scaling  Relationship 


An  empirical  expression  for  equilibrium  scour  depth  as  a  function  of  sediment  size  and 
maximum  discharge  per  unit  width  was  presented  in  CETN-IV-18  (Hughes  1999).  This 
relationship  was  formulated  as  a  balance  between  the  shear  stress  at  the  bottom  and  the  critical 
shear  stress  of  the  noncohesive  sediment.  An  unknown  coefficient  in  the  theoretical  formulation 
was  found  using  maximum  discharge  measurements  from  two  inlets.  The  resulting  formula  for 
equilibrium  scour  depth  is  restricted  to  those  portions  of  the  inlet  channel  where  currents  alone 
are  responsible  for  movement  of  sediment  by  bed  load.  The  equilibrium  scour  depth  relationship 
is  given  in  terms  of  discharge  per  unit  length  as: 


he  = 


0.234  4° 


8/9 


4/9  d'p 


(10) 


or  in  terms  of  velocity  as: 

V  =  5.1 24  [g  (Ss  -  l)]l/2  (11) 

where 


he  =  water  depth  at  equilibrium 
qe  =  maximum  equilibrium  discharge  per  unit  width 
g  =  gravitational  acceleration 
Ss  =  sediment  specific  gravity  [=ps  / pw  ] 
ps  =  mass  density  of  sediment 
pw  =  mass  density  of  water 
de  =  sand  median  grain  size 
V  =  depth-averaged  velocity 

Similitude  of  equilibrium  scour  depths  is  achieved  by  assuring  that  the  relationships  given  by 
Equations  10  and  11  have  the  same  effect  in  the  movable-bed  model  as  in  the  real  world.  The 
necessary  scaling  relationship  is  found  by  first  rearranging  Equation  11  into  the  form  of  a 
nondimensional  number,  i.e.: 


V 


[*(*.-!)] 


1,2  d, \K  A™ 


=  5.124 


(12) 


Similitude  requires  maintaining  the  same  value  of  the  nondimensional  number  in  both  the  model 
and  prototype,  or: 


\ 

V 

f  \ 

V 

{i8{sa-ir^ 

D 

oo 

oo 

cn  ^ 

<N 

H 

1 

>5 

Go 

(13) 
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where  the  subscripts  p  and  m  refer  to  prototype  and  model,  respectively. 

Equation  13  can  be  rearranged  into  ratios  of  prototype  values  of  a  variable  over  the 
corresponding  model  value.  These  are  termed  “scale  factors,”  and  they  are  denoted  by  an  upper 
case  N.  For  example,  the  scale  factor  for  depth- averaged  velocity  is: 


Replacing  all  the  prototype-to-model  ratios  with  scale  factors,  and  recognizing  that  gravity  will 
be  the  same  in  the  model  as  in  the  prototype,  results  in  the  following  scaling  relationship: 


Nv 


1 1/2 


=  1 


N, 


&"1) 


N**  NYS 


(14) 


or 


Nv  = 


n  1/2 


N, 


(^-1) 


n3/8  yyl/8 


(15) 


with  Nl  =  (he) p  /(he) m  representing  the  length  scale  of  the  model.  If  the  model  sediment  has  the 
same  density  as  the  prototype  sediment,  then  N(ss-i)  =  E  resulting  in  a  simplified  scaling 
relationship  of: 

Ny  =  N™  nY*  (16) 
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